Quasiperiodic localized oscillating solutions in the discrete nonlinear Schrodinger 

equation with alternating on-site potential 
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We present what we believe to be the first known example of an exact quasiperiodic localized stable 
solution with spatially symmetric large-amplitude oscillations in a non-integrable Hamiltonian lattice 
model. The model is a one-dimensional discrete nonlinear Schrodinger equation with alternating 
on-site energies, modelling e.g. an array of optical waveguides with alternating widths. The solution 
bifurcates from a stationary discrete gap soliton, and in a regime of large oscillations its intensity 
oscillates periodically between having one peak at the central site, and two symmetric peaks at the 
neighboring sites with a dip in the middle. Such solutions, termed 'pulsons', are found to exist 
in continuous families ranging arbitrarily close both to the anticontinuous and continuous limits. 
Furthermore, it is shown that they may be linearly stable also in a regime of large oscillations. 

PACS numbers: 42.65.Wi, 63.20.Pw, 63.20.Ry, 05.45.Yv 



The discrete nonlinear Schrodinger (DNLS) equation 
is one of the most studied examples of a non-integrable 
Hamiltonian lattice model. It is of great interest as well 
from a general nonlinear dynamics point of view, where 
it provides a particularly simple system to analyze fun- 
damental phenomena arising from competition of nonlin- 
earity and discreteness such as energy localization, wave 
instabilities etc., as from a more applied viewpoint de- 
scribing e.g. arrays of nonlinear optical waveguides 1] 
or Bose-Einstein condensates in external periodic poten- 
tials "^l . For recent reviews of the history, properties and 
applications of DNLS-like models, see Refs. 3, 4]. 

Recent attention has been given to DNLS-like models 
having, in addition to the fundamental periodicity given 
by the lattice constant, also a super lattice modulation 
creating a second period. In particular, the DNLS equa- 
tion with an additional term corresponding to alternating 
on-site energies has been proposed to model an optical 
waveguide array where the individual waveguides have 
alternating widths . Creating a two- branch linear dis- 
persion relation, the superperiodicity thus provides a pos- 
sibility for existence of a new type of nonlinear localized 
modes, discrete gap solitons (or discrete gap breathers), 
with frequencies in the gap between the upper and lower 
branches. These modes appear as stationary solutions to 
the DNLS equation (i.e., with constant intensity and a 
purely harmonic time-dependence which can be removed 
by transforming into a rotating frame), and their proper- 
ties were recently analyzed in detail in Ref. |^ , to which 
we also refer for further references on the topic. Very re- 
cently, they have also been experimentally observed 0- 

However, it is known that as the DNLS equation in 
one aspect is non-generic among nonlinear Hamiltonian 
lattice models, with a second conserved quantity being 
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the total norm of the excitation, there exist also local- 
ized quasiperiodic solutions which may have two (gen- 
erally) incommensurate frequencies. Here the first fre- 
quency corresponds to harmonic oscillations at constant 
intensity as for stationary solutions, while the second fre- 
quency corresponds to time-periodic oscillations of the 
intensity in the frame rotating with the first frequency. 
The existence of such solutions was proposed in Ref. 
and later explicit examples were constructed and ana- 
lyzed by continuation of multi-site breathers from the 
^anticontinuous^ limit of zero intersite coupling (note 
that similar ideas were used already in Ref. for finite 
systems) . A rigorous approach to the connection between 
existence of quasiperiodic solutions and additional con- 
served quantities was given in R ef. Illl| . A slightly differ- 
ent approach was taken in Ref. (see also 3]), where 
the existence of exact quasiperiodic solutions bifurcating 
from localized eigenmodes of the linearized equations of 
motion around some particular stationary solutions was 
shown. In all cases, to guarantee localization it is neces- 
sary that a non-resonance condition is fulfilled, so that no 
higher harmonics enter the continuous linear spectrum. 

However, all known examples of stable localized exact 
quasiperiodic solutions to the ordinary DNLS equation 
could be considered as rather special, since (i) they only 
exist in bounded parameter intervals at weak intcr-site 
coupling, and (ii) the intensity oscillations are typically 
quite small compared to their average values. Further- 
more, to the best of our knowledge no explicit example of 
a stable quasiperiodic breather with spatial symmetry has 
been given (although one of the modes presented in Ref. 
[l^ possibly could yield such a solution). On the other 
hand, it was shown already in Ref. 13] that for the two- 
dimensional DNLS model a state with large-amplitude 
symmetric intensity oscillations, with the intensity max- 
imum periodically oscillating between the central site 
and four surrounding sites, was created in an interme- 
diate stage in the process of 'quasicollapse' of a broad 
excitation to a highly localized breather. This entity. 
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termed ^pulson\ typically disappeared after 3-4 oscilla- 
tions, transforming into an on-site localized stationary 
breather with slowly decaying small-amplitude internal 
mode oscillations (the decay of which can be calculated 
similarly as for the one-dimensional case in Ref. [T^). 

It is the purpose of the present Brief Report to pro- 
vide the first example of an exact stable solution with 
such pulson properties, which we have found to exist 
in the one-dimensional binary modulated DNLS equa- 
tion. In fact, we remarked already in Ref. @ that sim- 
ilar arguments as in Ref. |l2l | prove the existence of ex- 
act quasiperiodic solutions bifurcating from certain in- 
ternal modes of the stationary gap breathers; here the 
crucial property (which is not fulfilled e.g. for single-site 
breathers in the nonmodulated DNLS model T^l) is, that 
these internal modes have frequencies above the continu- 
ous spectrum, and consequently also all higher harmonics 
will lie outside the continuous spectrum. Referring to the 
notation in Fig. 3 of Ref. we here concentrate on the 
spatially symmetric mode denoted as 'PI', which may 
exist arbitrarily close to as well the anticontinuous as the 
continuous limit, and show that the families of quasiperi- 
odic solutions bifurcating from this mode indeed exhibit 
pulson properties as the oscillation amplitude increases. 

Using for convenience a slightly different notation than 
in Ref. we consider the DNLS equation in the form: 

»^«-K)(-l)"V'n + C(^n+l+V'n-l) + IV'n|Vn =0, (1) 

with the two conserved quantities Hamiltonian H — 

En [Vo(-l)"|V'nP -C(V„C+1 +CV'n+l) " , 

and norm M = En IV'nP- Writing a stationary solution 
as V'n (t) = e*'^* with time-independent , the lin- 
ear dispersion relation becomes A = i + ^C^ cos^ q 
with gap A £ [— I Vol, I Vol]. Assuming Vo > 0, a discrete 
gap soliton then bifurcates from the lower gap edge 
A = — Vb, and corresponds in the limit C ^ to a single 
excited even site with intensity iV'nP = A + Vq- 

Exact quasiperiodic solutions with two independent 
frequencies, loq and Ub, can then be found numerically to 
computer precision by using essentially the same method 
as in Ref. 0. First, we transform into a frame rotating 
with the frequency c^o, '0n(i) = </'n(i)e"^°', yielding 

i<j)n~[uJ0 + Vo(-l)"] (t>n + C{(t)n+l+(t>n-l) + \4'nf4>n = 0, 

(2) 

whereupon standard Newton-type schemes are used to 
find time-periodic solutions to Eq. ((SJ fulfilling (j)n{t + 
27r/ti)b) — 0n(i), where is the second frequency. In ad- 
dition, linear stability is simply investigated numerically 
by standard Floquet analysis as described in Ref. P| . 

When a stationary gap breather of frequency A has an 
internal linear eigenmode of frequency lui (in the frame 
rotating with frequency A), the family of quasiperiodic 
solutions bifurcating from it at this point has ujb = loi 
and Wo = A — w;. Thus, by using as initial trial solu- 
tion the stationary gap breather perturbed with the rele- 
vant linear eigenmode, the two-parameter family of solu- 
tions (for given Vq and C) can then be followed through 



parameter space by standard continuation techniques. 
(Note that, in contrast to Ref. 9], once the stationary 
gap breather is known we here obtain directly quasiperi- 
odic solutions at finite C, without invoking the anticon- 
tinuous limit C — > 0.) For a small enough continua- 
tion step (< 10~^), we typically with standard double 
precision fortran routines for numerical integration and 
matrix inversion obtain convergence to solutions with 
\\(t)n{t + 2Tr/L0b) — 4'n{t)\\ < 10""^^ in five or less Newton 
iterations. In fact, since norm is a conserved quantity, it 
is generally more instructive to use TV and u>b as indepen- 
dent parameters, and practically continuations versus iOb 
at constant norm may be performed by varying ujq to keep 
TV constant. It is also convenient to work with rescaled 
quantities amounting to setting C = 1 in Eqs. ||2Jl; 
consequently the three relevant independent parameters 
for the two- frequency solutions can be chosen as Vq/C, 
TV/C (or, alternatively, ujq/C), and uJb/C. 

Fig. n illustrates the dynamics of a typical exact ex- 
ponentially localized pulson solution, obtained by con- 
tinuation (versus LOb/C at constant coq/C and Vq/C) of 
the linear 'PI' mode of a stationary gap breather with 
A = — 0.5Vo (i.e., with stationary frequency in the middle 
of the lower half of the gap) . The intensity oscillates pe- 
riodically, from having a minimum at the central (even) 
site and two symmetric maxima at the neighboring (odd) 
sites at t = 0, to having one single maximum at the cen- 
tral site a.t t ~ Tr/ujb- The solution is linearly stable, as is 
seen from the numerically calculated Floquet eigenvalues 
(Fig. n^b)), which are all on the unit circle. The main 
effect of a small but non-negligible perturbation is illus- 
trated by Fig.n(c): the solution almost perfectly retains 
its pulson character over very large time scales but with 
a slight shift in frequency (a consequence of the four- 
fold degenerate Floquet eigenvalue at -1-1 corresponding 
to drifts of the two arbitrary phases 0), yielding visible 
the pulson oscillations also in the stroboscopic plot. 

The results of more systematic investigations of the 
properties of such families of solutions are summarized 
in Figs. I2I4I For convenience, we now choose a 'moder- 
ate' constant value of the norm, TV/C = 2, and discuss 
the behaviour for various values of Vq/C. In particular, 
this value is large enough for the stationary gap breather 
frequency A to be sufficiently far away from the lower gap 
boundary so to always have a localized 'PI' mode, but 
small enough for A to be in the lower half of the gap, thus 
avoiding additional complications such as oscillatory in- 
stabilities and nonmonotonous continuation which rnay 
appear for stationary breathers in the upper gap half '(| . 

Fig. 12 for Vq/C — 2, shows the typical behaviour 
for larger values of Vq/C. The continuation of the two- 
frequency solution is monotonous, and ends in a bifurca- 
tion with another (unstable) stationary solution with fre- 
quency A — 2uJb + oJo- With the notation of Ref. j^, this 
solution is denoted as {O t (OO)}^, meaning that, in the 
anticontinuous limit C — > 0, it corresponds to symmetric 
in-phase oscillations only for the two (odd) sites ng ± 1, 
with frequency above the upper linear band. The pulson 
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FIG. 1: (a) The first three oscillations of an exact localized pulson solution to Eq. Q with Vo = 1, C = 1. The dynamical 
parameters are cub ~ 3.18, tJo = —3.722, TV « 1.90867. Only the central part is shown; the solution decays exponentially outside 
this region, (b) Linear stability eigenvalues for the solution in (a), (c) Snapshots, at integer multiples of the original period 
27r/tJi,, from long-time integration using the solution in (a), with a random perturbation ~ 10~^, as initial condition. 
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FIG. 2: Maximum and minimum amplitudes for the central 
(even) site no and its neighboring (odd) sites for continua- 
tion at constant norm A/" = 2 of the family of two-frequency 
solutions bifurcating from a stationary gap breather with 
A ~ —0.905 at ut ~ 4.186 (leftmost part of the figure). Here 
Vo — 2 and C — 1. The solution is unstable for cut > 4.47. 



character of the two- frequency solution appears when the 
absolute value of the minimum (negative) value of 4'no±i 
(lower part of dashed curve in Fig. O exceeds the mini- 
mum value of 4>nQ (lower part of solid curve), which for 
the case in Fig. |21 happens for ojb > 4.43. However, a sta- 
ble pulson solution appears in this case only in a small 
frequency interval, since for tob ^ 4.47 the solution is 
unstable through a symmetry-breaking instability (with 
real Floquet eigenvalue). For larger values oiVo/C, when 
the solution becomes more discrete, the solution becomes 
unstable before it attains its pulson character, and thus 
in the strongly discrete case the stable two-frequency so- 
lutions only correspond to relatively small oscillations of 
the central-site intensity for the stationary gap breather. 

For smaller values of Vq/C, the continuation versus ujb 
of the two- frequency solutions at fixed norm Af/C — 2 
still starts and ends in the same stationary solutions 
as before, but becomes nonmonotonous as shown for 
Vq/C = 1 in Fig. O (for clarity only minimum amplitudes 
(^„(0) are shown). The solution now attains its pulson 
character on the intermediate branch, where it is always 
stable, and becomes unstable only on the lower branch. 




FIG. 3: Same as Fig. |5|but for Vo = 1 and with only min- 
imum amplitudes. The two-frequency solution family bifur- 
cates from a stationary gap breather with A « —0.558 at 
uJb ~ 3.194, and is unstable for cut > 3.21 on the lower branch. 



For even smaller values of Vq/C, the stationary gap 
breathers broaden and assume a more 'continuous-like' 
shape, which in particular implies that the odd sites with 
largest l^i^^l will change from uq ± 1 to uq ± 3, and 
further to no ± 5, etc. |^. For Vq/C = 0.5 (Fig.©, the 
largest 'odd' \4'n'^\ for the stationary breather is at no ± 
5; however already for rather small intensity oscillations 
{ujb > 2.5856 on the upper branch of Fig. E^a)) of the 
corresponding 'PI' family of two- frequency solutions, the 
maximum 'odd' intensity is again at no ± 1. The solution 
now attains its pulson character already on the upper 
branch (ub ^ 2.5926), and the pulson remains stable until 
it reaches ujb ~ 2.5862 on the middle branch (Fig.^b)). 

Another very interesting effect regarding stability is 
seen m Fig. IHb). It is known 6] that as Vq/C de- 
creases, there are certain intervals of 'inversion of sta- 
bility' for stationary gap breathers, where the symmet- 
ric stationary breather discussed here becomes unstable 
(through a translational 'depinning' instability), and in- 
stead the otherwise unstable antisymmetric gap breather 
gains stability. Vq/C = 0.5 belongs to such an interval 
for Af/C = 2, and thus the stationary gap breather from 
which the two- frequency family in Fig. ^bifurcates is it- 
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FIG. 4: (a) Same as Fig. ^ but for Vo = 0.5. The two- 
frequency solution family bifurcates from the stationary gap 
breather with A ~ —0.3147 at cob ~ 2.5837, and (b) becomes 
unstable for uji, ~ 2.5862 on the middle branch. Magnification 
in inset of (b) illustrates how the weakly unstable stationary 
gap breather gets stabilized by the intensity oscillations. In 
(b), only the largest real part of the eigenvalue pair corre- 
sponding to the symmetry-breaking localized mode is shown. 

self unstable. However, as seen from the inset in Fig. 
^b), already at uib ~ 2.5856 the unstable eigenvalue re- 
turns to the unit circle, and thus even rather small sym- 
metric intensity oscillations may stabilize the stationary 



solution with respect to its antisymmetric instability. We 
are not aware of any reported similar scenario. 

To conclude, we have presented explicit examples of 
exact stable quasiperiodic pulson solutions with large- 
amplitude intensity oscillations, in the binary modulated 
DNLS equation describing e.g. coupled waveguides of al- 
ternating widths. Thus, there should be good possibili- 
ties for experimental observation of such states in this 
context. Although we here focused on one particular 
family of symmetric solutions bifurcating from the sym- 
metric stationary gap breather, the analysis illustrated 
by Fig. 3 in Ref. Q suggests that solutions with similar 
properties may bifurcate from other, symmetric or an- 
tisymmetric, internal modes with frequencies above the 
continuous spectrum ('P2', 'P3', 'P4', 'P5', etc.), ex- 
isting for the symmetric as well as the antisymmetric 
gap breathers. It seems likely, that such solutions also 
should exist for other types of multicomponcnt lattices 
with (at least) two conserved quantities; one interest- 
ing candidate being second- harmonic-generating lattices 
of the type considered e.g. in Ref. where also the 
fundamental discrete soliton was found to exhibit an in- 
ternal mode with frequency above the linear spectrum. 
A challenge for future research is to obtain analytical ex- 
pressions for the pulson solutions in the continuum limit. 
This is a nontrivial task, since although the relevant lin- 
ear eigenmodes of the stationary gap breather apparently 
persist arbitrarily close to the continuum limit (cf. Fig. 
3 of Ref. '^), they cannot exist in the standard contin- 
uous two-field model for gap solitons (which is the same 
for DNLS as for diatomic Klein-Gordon lattices, see e.g. 
Refs. [H 113), since the continuous spectrum of such 
models extends to infinity leaving no room for localized 
modes above it. A more sophisticated continuous approx- 
imation would therefore appear necessary to this end. 
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